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Wo study the decays B" D*°-k~ and B~ D*°K' , where the D*° decays into _D°7r°, with 
the D^ reconstructed in the CP-even (CP-I-) eigenstates K^K^ and 7r~7r+ and in the (non-CP) 
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channels K tt"*", K n'^n'^n , and K tv'^tt'^. Using a sample of about 123 million BB pairs, we 
measure the ratios of decay rates 



= 0-0813 ± 0.0040(stat)+°:--( 



and provide the first measurements of 



Rhp+ = = 0.086 ± 0.021(stat) ± 0.007(syst), 

B{B D^'pj^n ) 

and of the CP asymmetry 

A* = ^(B-^Dh%+K-) - B{B+^D*a%^K+) _ _q ^ q 23fstat)+» °^fsvst) 
AcP+ - B{B-^D'J>,+K-) + BiB+-.D'J>,+K+) " "'^^ ^ 0.23(stat)_o.o,(syst). 

PACS numbers: 14.40.Nd, 13.25.Hw 



The decays ^D^*^"K^*^~' will play an important 
role in our understanding of CP violation, as they can be 
used to constrain the angle 7 = aTg{-VudV*,,/VcdV*ij) 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix in 
a theoretically clean way by exploiting the interfer- 
ence between the b — > cus and b —>■ ucs decay am- 
plitudes In the Standard Model, neglecting D°I)0 
mixing, Rcp±/Rnon-CP — 1 + ''^ ± 2r cos (5 cos 7, where 
CP + (— ) indicates CP-even (odd) modes, 

p* _ _^ non-CP/CP± ' N 

^non-CP/ CP± = RCR- ^ r)*0 " T ' l^'' 

^non-CP/CP±^ > 

r is the absolute value of the ratio of the color suppressed 
B+ D*°K+ and color allowed B' D*°K- ampli- 
tudes (r ~ 0.1-0.3), and S is the strong phase difference 
between those amplitudes. The decays B^ D*^tt^ 
provide a convenient normalization term since many sys- 
tematic uncertainties are common to the two, while the 
interference effects should be highly suppressed for the 
D*^TT~ , when compared to the ones for the D*^K~ final 
states. Furthermore, defining the direct CP asymmetry 

B{B-^D-^%^K-)-B{B+^D*^%^K+) 
B{B-^D*^'p^K-)+B{B+^D*^p^K+y ^ ' 

we have: A^pi = ±2r sinJ sin7/(l -I- ± 2r cos (5cos7). 
The unknowns (5, r, and 7 can be constrained by measur- 
ing R*^on-CP^ ^cp±' ^cp±- The Belle Collaboration 
has reported RZ„^_cp = 0.078 ± 0.019 ± 0.009 using 10.1 
fb"^ of data 0. 

We present the measurement of -Rnon-CP' Rcp+ ^^"^ 
^CP+ performed using 113 fb~^ of data taken at the 
r(45') resonance by the BABAR detector with the PEP-II 
asymmetric B factory. An additional 12 fb~^ of data 
taken at a center-of-mass (CM) energy 40 MeV below the 
T{4:S) mass was used for background studies. The BABAR 
detector is described in detail elsewhere 0. Tracking of 
charged particles is provided by a five-layer silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH). 



The particle identification exploits ionization energy loss 
in the DCH and SVT, and Cherenkov photons detected 
in a ring-imaging detector (DIRC). An electromagnetic 
calorimeter (EMC), comprising 6580 thallium-doped Csl 
crystals, is used to identify electrons and photons. These 
systems are mounted inside a 1.5-T solenoidal super- 
conducting magnet. Finally, the instrumented flux re- 
turn (IFR) of the magnet allows discrimination of muons 
from other particles. We use the GEANT4 Monte Carlo 
(MC) j^J program to simulate the response of the de- 
tector, taking into account the varying accelerator and 
detector conditions. 

We reconstruct B^ — » D*'^h^ candidates, where the 
prompt track h~ is a kaon or a pion. Z)*° candidates are 
reconstructed from — > D'^tt'^ decays and mesons 
from their decays to K^ir^, X^tt^tt+tt", i^T^Tr+Tr*^, 
7r^7r+, and K^K^. The first three modes are referred to 
as "non-CP modes" , the last two as "CP modes" . Refer- 
ence to the charge-conjugate decays is implied here and 
throughout the text, unless otherwise stated. 

Charged tracks used in the reconstruction of D and 
B meson candidates must have a distance of closest ap- 
proach to the interaction point less than 1.5 cm in the 
transverse plane and less than 10 cm along the beam axis. 
Charged tracks from the 1?° — > 7r^7r+ decay must also 
have transverse momenta greater than 0.1 GeV/c and to- 
tal momenta in the CM frame greater than 0.25 GeV/c. 
Kaon and pion candidates from all D° decays must pass 
particle identification (PID) selection criteria, based on 
a neural-network algorithm which uses measurements of 
di?/dx in the DCH and the SVT, and Cherenkov photons 
in the DIRC. 

For the prompt track to be identified as a pion or a 
kaon, we require that its Cherenkov angle (Oc) be recon- 
structed with at least five photons. To suppress misrecon- 
structed tracks while maintaining high efficiency, events 
with prompt tracks with 9c more than 2 standard devi- 
ations (s.d.) away from the expected values for both the 
kaon and pion hypotheses are discarded; this selection 
rejects most protons as well. The track is also discarded 
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if it is identified with high probability as an electron or 
a muon. 

Neutral pions are reconstructed by combining pairs of 
photons with energy deposits larger than 30 MeV in the 
calorimeter that are not matched to charged tracks. The 
77 invariant mass is required to be in the range 122- 
146 MeV/c^. The mass resolution for neutral pions is 
typically 6-7 MeV/c^. The minimum total laboratory 
energy required for the 77 combinations is set to 200 MeV 
for 7r° candidates from D° mesons. Only 7r° candidates 
with CM momenta in the range 70-450 MeV/ c (denoted 
as soft pions, tTj) are used to reconstruct the D*^. 

The £)° mass resolution is 11 MeV/c^ for the — > 
K^tt^tt'^ mode and about 7 MeV/c^ for all other modes. 
A mass-constrained fit is applied to the D candidate. 
The resolution of the difference between the masses of the 
D*" and the daughter D° candidates (AM) is typically 
in the range 0.8-1.0 MeV/c^, depending on the Z)° decay 
mode. A combined cut on the measured D'^ and soft-pion 
invariant masses and on AM is also applied by means of 
a defined as: 





2 

+ 
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AM- AM 













2 
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where the mean values {rfijjo , rfiTr, , AM) and the res- 
olutions {om^Q, <^m„^, Cam) are measured in the data. 
Correlations between the observables used in the i^i 
Eq. (PJ are negligible. Events with > 9 are rejected. 

B meson candidates are reconstructed by combining a 
D*^ candidate with a high-momentum charged track. For 
the non-CP modes, the charge of the prompt track h must 
match that of the kaon from the meson decay. Two 
quantities are used to discriminate between signal and 
background: the beam-energy-substituted mass tties = 
^y{E*'^/2 + Pi ■ tpbY /E'f - pg and the energy difference 
iS.E = Eg — E* /2, where the subscripts i and B refer to 
the initial e^e~ system and the B candidate respectively, 
and the asterisk denotes the CM frame. 

The mES distribution for B~ D*^h~ signal can be 
described by a Gaussian function centered at the B mass 
and does not depend on the nature of the prompt track. 
Its resolution, about 2.6 MeV/c^, is dominated by the 
uncertainty of the beam energy and is slightly dependent 
on the i?" decay mode. The observable AE does depend 
on the mass assigned to the tracks forming the B candi- 
date, and on the momentum resolution. We calculate 
AE with the kaon hypothesis for the prompt track and 
indicate this quantity with AEk- For B~ D*°K^ 
events AEk is described approximately by a Gaussian 
centered at zero and with resolution 17-18 MeV, whereas 
for B~ D*'^TT~ events AEk is shifted positively by 
about 50 MeV. B candidates with tues in the range 5.2- 
5.3 GeV/c^ and with AEk in the range (-100 to 130) 
MeV are selected. 

A large fraction of the background consists of contin- 



uum (non BB) events and a powerful set of selection cri- 
teria is needed to suppress it. The selection is chosen to 
maximize the expected significance of the results, based 
on MC studies. In the CM frame, this background typi- 
cally has two-jet structure, while BB events are isotropic. 
We define 9t as the angle between the thrust axes of 
the B candidate and of the remaining charged and neu- 
tral particles in the event, both evaluated in the CM 
frame, and signed so that the thrust axis component 
along the e~ beam direction be positive. The distribu- 
tion of I cos 9t I is strongly peaked near one for continuum 
events and is approximately uniform for BB events. For 
the non-CP modes, |cos0t| is required to be less than 
0.9 for the K~Tr'^ mode, and less than 0.85 for 

D'^ if ^TT+TT+TT^and D'^ K^ir^n^ modes for which 
the levels of the continuum background are higher. For 
the CP modes, cos 6t is required to be in the ranges (—0.9 
to 0.85) and (-0.85 to 0.8) for the £>" R- K+ and 
— > TT^TT"*" modes respectively. Other mode-dependent 
selection criteria are applied: for the _fC~7r+7r+7r~ 
and D'^ K^n^n^ {D^ 7r^7r+) modes we reject 
events with cosOtD < —0.9 (|cos0(£i| > 0.95), where 
OtD is the angle between the direction of the in the 
laboratory and the opposite of the direction of the 
{tt~ for the 7r~7r"'" mode) from the D'^ in the 

rest frame. Finally, to reduce combinatorial background 
in the K^tt^t:'^ final state, only those events that 

fall in the enhanced regions of the Dalitz plots, according 
to the results of the Fermilab E691 experiment J^], are 
selected. This last requirement alone rejects 80% of the 
background and accepts 69% of the signal, according to 
the MC simulation. 

Multiple candidates are found in about 10-12% of the 
selected events with two- and four-body D'^ decays and in 
17% of the events with D° K-Tr+n° decays. The best 
candidate in each event is selected based on the Pre- 
viously defined. The number of candidates constructed 
with the same D*", but different prompt track, is negli- 
gible; in this rare case the best one in the event is ran- 
domly chosen. The reconstruction efhciencies, based on 
MC simulation, are reported in Tabled 

According to the simulation, the main contributions 
to the BB background for B^ D*'^h~ events origi- 
nate from the decays B^ ^ D(*)°p^ and B° D*~h+ . 
B^ D*'-\-^ D'^j)h^ events are also considered back- 
ground as their CP modes have CP eigenvalues opposite 
to the ones of the B~ — > D*'^ h~ signal 

For each decay mode, an unbinned maximum- 
likelihood (ML) fit is used to extract yields from the data 
for six candidate types: signal, continuum background, 
and BB background, for the kaon and pion choices for 
the mass hypothesis of the prompt track in the candidate 
decays B' D*°h- . 

Three quantities from each selected candidate are used 
as input to the fit: AEk, wesi and the 9c of the prompt 
track. The distributions of AEk and mEs for the six 
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candidate types are parametrized to build the probability 
density functions (PDFs) that are used in the fit. 

Correlations between the jties and IS.Ek variables for 
signal events are about —5% according to the simula- 
tion. To account for these, we use signal MC events to 
parametrize the signal PDFs with a method based on 
kernel estimation Q, which allows the description of a 
two-dimensional PDF. The shapes of MC and data dis- 
tributions of these observables are in good agreement, 
according to comparisons performed with pure samples 
of B~ D*^TT~ events, obtained with very tight particle 
identification and kinematic selection. To the extent that 
we find differences in the data and MC distributions, we 
adjust the shapes of the PDFs to conform to the data. 
Systematic uncertainties due to limited statistics associ- 
ated with this procedure are included in the final results. 

We obtain the PDFs for the mES distribution for 
continuum background from off-resonance data, apply- 
ing the standard selection criteria. The toes distri- 
butions are parametrized with a threshold function 
defined as /(toes) oc y\/l — y"^ exp [—^(1 — y^)], where 
y = toes/wq and toq is the mean energy of the beams 
in the CM frame. The PDFs for the AEk distribu- 
tions for background candidates from the continuum are 
well parametrized with exponential functions whose pa- 
rameters are determined by fitting the AEk distribu- 
tions of the selected B~ — > D*^h~ sample in the off- 
resonance data. Both the toes and the AEk PDFs for 
the continuum background are taken to be the same for 
B- D*°Tr- and B- D*°K- decays. The shapes 
of MC and data distributions of toes and AEk obtained 
with looser selection criteria to increase the statistics, 
agree well for B^ D*'^tt^ and B^ D*^K^ decays, 
validating this assumption. For the CP modes very few 
off-resonance events pass the selection criteria, hence we 
use the PDFs determined for the — * A'^7r+ mode. 
This is justified by a separate comparison of the CP 
modes with the flavor-definite modes in data and MC 
samples obtained with looser selection criteria. 

The correlation between toes and AEk for the BB 
background is taken into account with a two-dimensional 
PDF determined from simulated events, in a similar way 
to that used for the signal. 

We obtain PDFs for the particle identification deter- 
mination for the prompt track from the distributions, in 
bins of momentum and polar angle, of the difference be- 
tween the reconstructed and expected 9c of kaons and 
pions from decays in a control sample that exploits 
the decay chain D*+ D°tt+, K-tt+ to identify 

the tracks kinematically. 

Initial PDFs are parametrized for each candidate type 
as detailed above. With these we then fit pure samples 
of simulated signal events and of background from off- 
resonance real and MC data. With the yields from these 
fits we establish an efficiency matrix accounting for small 
crossfeeds among the components. The corrections af- 



fecting the signal yields are typically of order 1%. The 
fractional systematic uncertainties for the signal yields 
associated with these corrections are in the range 0.1- 
6.0% depending on the decay mode. 

The likelihood C for the selected sample is given by the 
product of the final PDFs for each individual candidate 
and a Poisson factor: 

nE^^.-("^BS.,Ai^^,^cJ (4) 

i=lj=l 

where N is the total number of events, Nj are the yields 
for each of the previously defined six candidate types, 
and N' = VjiniESi, AEk,, 9ci) is the proba- 

bility to measure the particular set of physical quanti- 
ties {mESi:AEKi,Oci) in the z'^ event for a candidate of 
type j. The Poisson factor is the probability of observ- 
ing N total events when N' are expected. The quantity 
C is maximized with respect to the six yields using the 
MINUIT program 0. The fit has also been performed 
on luminosity-weighted MC and high statistics toy MC 
events and it has been found to be unbiased. 

The results of the fit are reported in detail in Table 
These yields are used to determine the CP asymmetry 
parameters. We measure: 

Kon-cp = 0.0813 ±0.0040(stat)+°;°°^?(syst), 

Rcp+ = 0.086 ±0.021(stat)±0.007(syst), 

Rhp+/Kon-CP = 1.06±0.26(stat)+°:;°(syst), 

^CP+ = -0.10±0.23(stat)+o°4(syst). 

Figure m shows the distributions of AEk for the com- 
bined non-CP and CP modes before and after the en- 
hancement of the B D*^K component. The en- 
hancement is accomplished by requiring that the prompt 
track be consistent with the kaon hypothesis and that 
Toes > 5.27GeV/c^. The AEk projections of the fit re- 
sults are also shown. 

TABLE I: Results of tlie yields from the ML fit. For the 
CP modes the results of the fit separately for the B'^ and 
B~ samples are also quoted. Errors are statistical only. The 
efficiencies (e) based on MC simulation are also reported. 



i?" mode 


N{B D'^n) 


N{B D*°K) 


£(£)*%) (%) 


K'-!V + 


2639 ± 56 
3249 ± 68 
3071 ± 64 


226 ± 18 
247 ± 21 
242 ± 21 


17.5 ±0.2 
5.9 ±0.1 
9.7 ±0.1 


K-K+ 
K-K+ [B+] 
K-K+ [B \ 


258 ± 19 
123 ± 13 
134 ± 13 


23.4 ±5.6 
13.4 ±4.1 
9.9 ±3.7 


15.3 ±0.2 
15.6 ±0.3 
14.9 ±0.3 


■K-TT+ [B + ] 
■K--K+ [B \ 


124 ± 14 
75 ±11 
49 ±9 


6.3 ±4.6 
0.7 ±3.2 
5.3 ±3.5 


14.6 ±0.2 
14.5 ±0.3 
14.8 ±0.3 
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FIG. 1: Distributions of \Ek in the B -> sam- 
ple, for D° ~* K-'-!v+,K--!v+TT°,K-TT+Tv+-!v~ ((a), (b)) and 
D° K-K+,Tv-TT+ ((c), (d)), before ((a), (c)) and after 
((b), (d)) enhancing the B D*^K component by requiring 
that the prompt track be consistent with the kaon hypothesis 
and mES > 5.27 GeV/ The B" ^ D"-'' n signal contribu- 
tion on the right of each plot is shown as a dashed line, the 
B~ D*°K~ signal on the left as a dotted line, and the 
background as a dashed-dotted line. The total fit with all the 
contributions is shown with a thick solid line. 



The ratio of the decay rates for D*'^Tr^ and 

D*^K'' is separately calculated for the different 
decay channels and is computed with the signal yields 
estimated with the ML fit and listed in Table HI The re- 
sulting ratios are scaled by correction factors of a few 
percent, which are estimated with simulated data and 
which take into account small differences in the efficiency 
between B~ — » D*°K~ and B~ — > D*°Tr~ event selec- 
tions. The results are listed in Table Hll 

TABLE II: Measured ratios for different _D° decay modes. 
The first error is statistical, the second is systematic. 



B- D*°h- Mode 


B(B^D*°K)/B{B~>D''\) (%) 


D« ^ K--K+ 


8.93 ± o.72i;^';^* 

7.59 ± 0.65tg;^^ 
7.91 ± 0.72115 


Weighted Mean (non-CP) 


8.13 ±0.40+0 31 


D° K-K+ 
D° ^ 7r-7r+ 


9.4 ±2.3 ±0.6 
5.9±4.4l:l;4 


Weighted Mean (CP) 


8.6 ±2.1 ±0.7 



The sources of systematic uncertainties for the yields 
have been identified and their contributions (for the mea- 
surement of R*non-)Cp) reported in Table HTll Uncer- 
tainties of the signal parametrizations of AEk and toes 



TABLE III: Average systematic uncertainties for 



Systematic 


^-^iioii-CP /-^non-CP (^) 


AR*cp/R'cp(%) 


Source 


non-CP modes 


CP modes 


AEk (signal) 


+2.0 
-1.8 


+2.7 
-2.8 


A VP rtri\ 

L\tjK\qq) 


+0.3 
-0.6 


+0.9 
-2.5 


AEk{BB) 


+0.0 
-0.5 


+0.8 
-0.8 


mES (signal) 


+0.4 
-0.3 


+0.7 
-0.4 


niES (qq) 


+0.8 
-0.8 


+4.4 
-6.7 


mES (BB) 


+ 1.2 
-1.3 


+0.3 
-3.2 


PDF Crossfeeds 


+2.8 
-2.6 


+0.7 
-0.7 


PID PDF 


+3.0 
-1.8 


+4.0 
-1.4 


£ Correction 


+ 1.4 
-1.4 


+2.0 
-2.0 



arise from the assumed shapes of the PDFs and discrep- 
ancies between real and simulated data. All the param- 
eters of the AEk and toes PDFs have also been varied 
according to their one s.d. statistical uncertainties and 
signed variations in the yields are taken as systematic 
uncertainties. For the BB and continuum backgrounds, 
the systematic uncertainties due to the limited statistics 
of the MC and of the off-resonance data have been cal- 
culated varying the AEk and toes PDF parameters by 
their statistical uncertainties. There are several contribu- 
tions to the PfD systematic uncertainty for the prompt 
track: the uncertainty due to limited statistics is calcu- 
lated by varying each parameter of the PDF, in each bin 
in momentum and polar angle, by its uncertainty (keep- 
ing constant all other parameters in the same bin and 
all parameters in all the other bins) and summing all the 
contributions in quadrature; results obtained with alter- 
native PID PDFs, which account for different 6c residual 
shapes and for discrepancies between data and simula- 
tion, are also included as systematic uncertainties. The 
systematic uncertainties due to the fit crossfeeds have 
been evaluated. Finally, errors associated with the effi- 
ciency correction factor are also included. 

Many of the systematic uncertainties for the signal 
yields have similar effects on the B^ D*^K~ and 
B^ D*^TT^ events (they increase or decrease both 
fractions simultaneously), hence their effect is reduced 
in deriving the systematic uncertainty for the measure- 
ment of the ratios, when all correlations are taken into 
account. Overall, the main sources of systematic un- 
certainties for the measurement of both i?T ^„r, and 

( noil- J 

^CP+ ^'^'^ '^^'^ ^'^ characterization of the shapes of 
Toes and AEk for the signal, to the characterization of 
the Toes PDFs for the background, to the particle iden- 
tification, and to the uncertainty of the fit crossfeeds and 
of the efhciency correction factors. The systematic un- 
certainty for A*(jpj^ due to possible detector charge asym- 
metries is evaluated by measuring asymmetries analogous 
to those defined in Eq. |(2Jl, but for B~ — > D*^tt~ and 
B^ D*'^K^ events (the latter uniquely for the non- 
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CP modes), where CP violation is expected to be negli- 
gible. Results for all modes are then combined, taking 
correlations into account. The measured asymmetry is 
-0.008 ± 0.012(stat) ± O.OOl(syst). Though it is consis- 
tent with zero, it is also consistent with —0.020 at one 
s.d. level, hence we take the magnitude of this value 
as a further symmetric systematic uncertainty on A^p^ . 
When combining the results for the different modes, all 
systematic and statistical uncertainties are considered to 
be uncorrelated, except for the contributions of the PID 
PDF (common to all modes) and of the detector charge 
asymmetry in the measurement of ^Jy^ , which are con- 
sidered to be completely correlated. For the measure- 
ment of Rcp+I Rnon-CP Systematic uncertainties have 
been considered to be uncorrelated; this assumption is 
conservative, and has negligible effect on the final result, 
which is largely statistically limited. 

In conclusion, we have measured the ratio of the de- 
cay rates for B~ D*°K~ and B~ D*^tt^ pro- 
cesses, with non-CP eigenstates. This constitutes the 
most precise measurement for this channel. We have 
also performed the first measurement of the same ra- 
tio and of the CP asymmetry ^cp+ for D° mesons 
decaying to CP eigenstates. These results, together 
with measurements exploiting B~ D^K~, B~ 
D"K*- and ^ D*°K*- decays 0, constitute 
a first step towards measuring the angle 7. Further- 
more, assuming factorization and flavor-SU(3) symme- 
try, theoretical calculations (in the tree-level approxima- 
tion) predict: B{B- D*°K-)/B{B- D*°tt-) - 
{Vus/Vndy{fK/U f ~0.074, where fa and/^ are the me- 
son decay constants Our results accord with these 
predictions. 
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